Unstable growth of curvature perturbation in non-singular bouncing cosmologies 
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We consider non-singular bouncing cosmologies, such as the new ekpyrotic model, in which the 
universe undergoes a slow contraction phase with equation of state w ^ 1, followed by a bounce 
that occurs at a finite scale factor when quantum gravity corrections are still negligible. Such a 
non-singular bounce requires a violation of the null energy condition in which w falls below —1 at 
some time before the bounce. In this paper, we show that a component of the adiabatic curvature 
perturbations, though decaying and negligible during the ekpyrotic phase, is exponentially amplified 
just before w approaches —1, enough to spoil the scale-invariant perturbation spectrum. We discuss 
how the problem may be avoided, for example, in singular bounces. 



o 



I 



> 

00 
(N 

l> 

O 
O 



X 



As an alternative to inflationary cosmology, the ekpy- 
rotic scenario [H |2] resolves the homogeneity, flatness 
and isotropy problems of standard cosmology by hav- 
ing the universe undergo a period of ultra-slow contrac- 
tion followed by a bounce to the standard expanding 
phase. Ultra-slow contraction [ekpyrotic phase) occurs 
when the dominant energy component has an equation of 
state It; 3> 1. During contraction, this component grows 
faster than all other contributions to the cosmic expan- 
sion, including the spatial curvature and anisotropy, and 
thereby drives the universe into an exponentially flat and 
isotropic state leading up to the bounce ^3j. 

As for the bounce, two possibilities have been dis- 
cussed. In "singular bounce" models, such as the cyclic 
model [4| , the universe contracts towards a "big crunch" 
until the scale factor a{t) is so small that quantum grav- 
ity effects become important. The presumption is that 
these quantum gravity effects introduce deviations from 
conventional general relativity and produce a bounce that 
preserves the smooth, flat conditions achieved during the 
ultra-slow contraction phase. 

The second possibility is a "non-singular bounce" , such 
as in the "new ekpyrotic model" [5], where the universe 
stops contraction and reverses to expansion at a finite 
value of a(t) where classical general relativity is still valid. 
A significant advantage of this scenario is that the entire 
cosmological history can be described by 4d effective field 
theory and classical general relativity, without invoking 
extra dimensions or quantum gravity effects. However, 
for the bounce to happen within classical general relativ- 
ity, w must fall below —1 for a sustained period, i.e. a 
violation of the null energy condition (NEC). The NEC 
violating energy component is commonly chosen to be a 
scalar field that undergoes ghost condensation [51 [Z]- 

In this paper, we show that the non-singular bounce 
creates a serious problem for cosmological perturba- 
tions. Observations of the cosmic microwave background 
(CMB) and large-scale structure have provided evidence 
for a scale-invariant perturbation spectrum. However, 
in the ekpyrotic model when the scale-invariant curva- 
ture perturbation is generated during or just after the 
ekpyrotic phase, a potentially dangerous component of 



adiabatic curvature perturbations is created at the same 
time. This mode has been previously ignored because, 
after exiting horizon when w > 1, its amplitude becomes 
exponentially suppressed on large length scales compared 
to the scale-invariant modes. In a singular bounce, this 
mode remains completely negligible because w > 1 all 
the way up to the bounce. However, for the non-singular 
bounce, the ekpyrotic phase must end and w must fall 
below —1. We show that, right before crossing w — —1, 
the adiabatic mode undergoes exponential amplification 
such that the scale-invariant spectrum is spoiled and per- 
turbation theory breaks down. 

For our study we take the new ekpyrotic model [5] 
as a generic example of non-singular bounces. In this 
model both the ekpyrotic phase and the ghost-condensate 
phase are driven by the same scalar field, since otherwise 
one encounters the problem that the energy of the ghost- 
condensate field is so exponentially diluted away after the 
ekpyrotic phase that it is too insignificant to cause the 
bounce. This framework can be described by an effective 
Lagrangian 
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{\R + P{X)~V{c^)), X^~\{dc^)\ (1) 



assuming a FRW background metric with signature ( — h 
-I— 1-), and reduced Planck units SttG = 1. The kinetic 
term P{X) is designed as in Fig. [l] where it is linear for 
large X, P{X) « X, but has a minimum at a low en- 
ergy scale Xc. The potential V{(l)) is sketched in Fig. [5] 
where, beginning from the right hand side, V is approx- 
imated by a negative exponential — Voe^v^/J"^ over a 
range between T4k-beg and T4k-cnd, then bottoms out and 
undergoes a steep rise. The universe evolves through the 
ekpyrotic, kinetic and bouncing phases, as indicated in 
the figures. 

During the ekpyrotic phase, the Friedmann equation 
has an attractor solution. 
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where t is negative and increasing towards zero, and we 
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FIG. 1: The kinetic term P(X) versus X = \{d(l>f for the 
ghost-condensate scalar field (j). As indicated by the dotted 
lines with arrows, during the ekpyrotic phase {X > Xc), the 
approximately canonical kinetic term P{X) ~ X grows expo- 
nentially by a factor of e^'^'"' ; After the ekpyrotic phase, X 
decreases by an even greater factor to reach X = Xc. This 
exponential decrease in X is directly related to the problem 
with non-singular bounces. 
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FIG. 2: The effective potential V{<j}) in the ekpyrotic model 
with a non-singular bounce. The evolution along the poten- 
tial during the ekpyrotic and bounce phases is from right to 
left: the ekpyrotic phase refers to the exponential decline from 
V^k-bcg to l/ek-cnd near the minimum of the potential; then 
there is a brief kinetic energy dominated phase as V quickly 
rises; and finally the non-singular bouncing phase occurs as 
V climbs sufficiently above zero. 



normalize flok-cnd = 1- This solution has a constant equa- 
tion of state, w — J — 1^1, if we choose p <^ 1. The 
potential energy and the nearly canonical kinetic energy 
satisfy the scaling relation 
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Therefore, there is an exponential suppression of any ini- 
tial curvature and anisotropy, determined by the ratio of 
the kinetic energy densities at the end and beginning of 
the ekpyrotic phase. 
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A brief kinetic energy dominated phase follows after 
the field reaches the bottom of the potential and rises 



towards V = 0. During this phase the equation of state 
w rapidly decreases according to 
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The first term in the brackets represents the ratio of the 
gradient force to the total energy, which is a large positive 
factor during the sharp rise in the potential. Therefore 
w < in this phase, since H is negative. 

In singular bounces [3], the potential approaches zero 
from below as </> — ^ — oo, so that w remains > 1 until 
the bounce. In contrast, for a non-singular bounce, the 
potential rises above zero as in Fig. [2] so that, as the 
field climbs, w decreases below 1 and eventually crosses 
— 1 some time before the bounce. This stage occurs in 
much less than one Hubble time, during which H and a 
are nearly constant. The kinetic term in the Lagrangian 
remains canonical for X ^ X^. 

The bouncing phase begins when the field climbs suffi- 
ciently far up the potential that X falls below X^ and 
the ghost condensate phase initiates. In this phase, 
H = —XP^x = ^5P(1 + w) becomes positive and w falls 
below —1, violating the NEC. The non-singular bounce 
occurs as contraction slows and ultimately halts at a fi- 
nite value of a. 

For consistency, a necessary condition is (Fig. IlJ 
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which keeps the kinetic term P{X) canonical throughout 
the ekpyrotic phase. The speed of sound is 
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when P{X) is linear; in the bouncing phase when w < 
— l,c^ becomes small and negative, but it does not cause 
instabilities if the universe bounces and exits the ghost- 
condensate phase within a Hubble time or so [SJ |H] . 

The serious problem appears when we consider per- 
turbations around the background evolution. Here we 
focus on the widely studied gauge-invariant variable, the 
comoving curvature perturbation TZ. Perturbations in 
other gauges will be discussed in [5] , where we show that 
results agree at the bounce, so that the initial conditions 
for the expanding phase are consistent in all gauges. Dif- 
ferent Fourier modes TZk of the comoving curvature per- 
turbation, labeled by the comoving wave number k, obey 
the equation of motion ^ 
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where z = ci\/ -^jj2, and prime denotes derivative with 

respect to conformal time r, ^ = o^- For small k, the 
equation is formally solved in expansions of k^ , 

, (0) ,-, f dr 
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where the leading order term is the sum of a constant 
term T?.™"**' and an integral term T?."* 
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The diniensionless constants Ci and C2 are numbers of 
order 0(1), as can be found by matching to Minkowski 
vacuum condition at the beginning of the ekpyrotic phase 
when the mode is deep inside the horizon. 

During the ekpyrotic phase, ^ = |(1 + w) = ^ 
is nearly constant and a ^ {—t)^ does not change sig- 
nificantly. Therefore, the integral term decreases as 
^fe"* '^ H-t) ast-)-0-. It is comparable to 7^^°"''' at 
horizon crossing (fc ^ aH), but becomes exponentially 
small by the end of the ekpyrotic phase. 
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where X^ is the kinetic energy at horizon crossing and 
Nk is the remaining number of c-foldings of the ekpy- 
rotic phase after k mode exits horizon. (Factors of p are 
neglected for the purpose of these estimates.) Thus 
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Unfortunately, 7?.™"''* has a blue spectrum {Pji (x 
fc3|7^fe|2 ex k^), and TZ^* is bluer still (P-r oc k^), in- 
consistent with the scale-invariant spectrum observed in 
CMB. 

Following |S] , we consider the entropic mechanism [10] 
for generating a scale-invariant spectrum. By introduc- 
ing two scalar fields ^i and 02 that undergo an ekpyrotic 
phase simultaneously, we have an extra degree of freedom 
that can source the curvature perturbation. Fluctuations 
in the fields are decomposed into an adiabatic mode da 
along their mean trajectory and an entropic mode 6s per- 
pendicular to it. Both fields obtain scale- invariant fluc- 
tuations during the ekpyrotic phase, but they source the 
curvature perturbation 7?. differently [T^, through the 
equation 
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The first term represents the adiabatic contribution TZ^°'' 
to the curvature perturbation, which is the same as in 
the single field case. The second term represents the 
entropic contribution TZ^'^' , where 6 is the angular di- 
rection of the trajectory in the {(pi, 412) plane. It is as- 
sumed that the trajectory is nearly straight except 
at the end of the ekpyrotic phase, where it undergoes a 
sharp bend and renders 9 temporarily non-zero. Accord- 



ing to Eq. (13), this causes the entropic perturbation 6s 



to convert almost instantaneously into a scale-invariant 



curvature perturbation 7?.'*' , which remains constant on 
superhorizon scales afterwards. 

Thus, the total curvature perturbation can be decom- 
posed as 

-^tot ^ j^{s) ^ ^{a) _ ^(s) ^ ^(,T,const) ^ 7^(.T,int)^ j^j^-j 

where as before the adiabatic contribution is divided into 
a constant term that is blue and an integral term that 
is bluer. For the modes that exited horizon during the 
ekpyrotic phase, the constant term T^lc'jconst) jg subdom- 
inant compared to the scale-invariant contribution TZ^''', 
and the integral term TZ^°''™^' is sub-subdominant, 



(15) 
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The critical stage occurs towards the end of the kinetic 
phase after the conversion. While 7^('') and 7^('^.™nst) 
remain constant, T?,'"^'™'' grows exponentially as w — >■ — 1. 



The growth can be seen by rewriting Eq. ( 10 1 as 
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Integrating w from lu 3> 1 at the end of the ekpyrotic 
phase to ui ~ — 1 just before the non-linearity in P{X) be- 
come significant, the dominant part of the integral comes 
from near w w — 1. Using Eq. ([5]) and neglecting the sec- 
ond term as w approaches —1, we obtain 
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where we have neglected some finite factors that are al- 
most constant during the rapid decrease in w, and used 
relation ^ since P{X) is still linear. The integral term 
has grown exponentially compared to its value at the end 



of the ekpyrotic phase from Eq. ( 12 ) 
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Hence, from Eq. ( 16 ), the ratio of the integral term to the 



scale-invariant entropic contribution can be expressed as 
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This ratio determines whether the adiabatic contribu- 
tion can catch up with the entropic contribution and 



dominate the curvature perturbation. There is a com- 
petition between two factors: (i) an exponential sup- 
pression from ( 16 1 that depends on how much X in- 
creases after horizon exit during the ekpyrotic phase, 
i.e. Xk/Xc\i-cnd ~ e~^^* as in Eq. (11); and (ii) an ex- 
ponential amplification from ( 19 ) that depends on how 
much X decreases in the kinetic phase, i.e. ^ok-ond/^c > 
-'^ok-ondZ-'^ek-beg = 6^^'°' according to Eq. Q. Therefore 
we find, in total. 
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The inequality boils down to the fact that the increase of 
the kinetic energy X in the ekpyrotic phase is less than 
the decrease in the kinetic phase, hence amplification 
wins: This ratio is much greater than unity for modes 
with Nf: < iVtot/2, which includes all modes within our 
observable horizon. That is, the integral adiabatic con- 
tribution, which has a doubly blue spectrum and was 
sub-subdominant at the end of the ekpyrotic phase, has 
grown to overwhelm the scale-invariant entropic contri- 
bution on relevant scales by the time the bouncing phase 
begins. 

Thus, the dominantly blue spectrum will be carried 
into the expanding phase, in contradiction to current 
observations. Moreover, as Eq. (21) dictates, for a 



wide range of k modes the perturbation amplitudes have 
grown non-linear, to the extent that perturbation theory 
may break down even before the bounce. 

As captured in Fig. [l] the problem with the non- 
singular bounce arises from the generic requirement that 
in the kinetic phase the kinetic energy density must de- 
crease by more than it has increased in the ekpyrotic 
phase, in order to trigger the ghost-condensate phase and 
violate the NEC. Hence, it is difficult to avoid this prob- 
lem in a ghost-condensate bouncing model with a canon- 
ical (linear P{X)) ekpyrotic phase that does not involve 
extreme fine-tuning. Although the hope had been that a 
non-singular bounce would make ekpyrotic models sim- 
pler than the singular case by not having to consider 
quantum gravity effects, we are forced to conclude that 
this scenario leads to the exponential growth of adiabatic 
perturbations that spoil the spectrum of density pertur- 
bations. 

As a possible remedy, the growth of curvature per- 
turbations may be moderated if the ekpyrotic phase is 
realized with some non-linear form of P{X,<j)), e.g. in 
[131 [T3]. This kind of non-linearity is typically associ- 



ated with quantum gravity effects, as occurs naturally in 
singular bounces. Hence, the simple non-singular bounce 
may fail for the reasons described in this paper, but it re- 
mains possible that a non-linear realization |14| , or a sin- 
gular bounce, can produce an observationally acceptable, 
nearly scale-invariant spectrum of density perturbations. 
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